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In this work we develop a numerical framework to investigate the renormalization of the non- 
Markovian dynamics of an open quantum system to which dynamical decoupling is applied. We 
utilize a non-Markovian master equation which is derived from the non-Markovian quantum trajec- 
tories formalism. It contains incoherent Markovian dynamics and coherent Schrodinger dynamics as 
its limiting cases and is capable of capture the transition between them. We have performed com- 
prehensive simulations for the cases in which the system is either driven by the Ornstein-Uhlenbeck 
noise or or is described by the spin-boson model. The renormalized dynamics under bang-bang 
control and continuous dynamical decoupling are simulated. Our results indicate that the renormal- 
ization of the non-Markovian dynamics depends crucially on the spectral density of the environment 
and the envelop of the decoupling pulses. The framework developed in this work hence provides an 
unified approach to investigate the efficiency of realistic decoupling pulses. This work also opens a 
way to further optimize the decoupling via pulse shaping. 



I. INTRODUCTION 

It is known that dynamical decoupling 0,0,0)01 is one 
of the essential control tools in fighting the decoherence 
of an open quantum system. In recent years many meth- 
ods have been developed to design the decoupling pulses. 
Those methods are based on the group structure or the 
geometric perspective of the system environment interac- 
tion '3', '4', '5', '3, T, 'sj . The effect of a general environment, 
however, is usually neglected. Dynamical decoupling can 
be divided into several categories. It can be determinis- 
tic and periodic'!','^, deterministic and concatenated[Tfil| . 
or random 8, 11]. The strength of the pulses can be un- 
bounded, which is usually termed bang-bang decoupling, 
or bounded 0, 0| ■ Bounded control uses more realistic 
control resources and is more tolerant against implemen- 
tation errors. Both unbounded and bounded control at- 
tain the desired decoupling in the ideal limit where the 
period of a single decoupling cycle approaches zero. In 
practice, however, the duration of the cycle is finite and 
there will always be residual decoherence. To the best of 
our knowledge it is difficult to calculate the exact resid- 
ual decoherence, especially when the effect of a general 
environment are taken into account. For a general config- 
uration the amount of residual decoherence might depend 
on the shape of the decoupling pulses as well as the spec- 
tral density of the environment. For example, it has been 
argued that the bang-bang decoupling might be more ef- 
ficient when the environment is characterized by the l/ f 
noise in stead of an Ohmic bosonic environment [13. IT^ . 
Recently the effect of a general environment is studied, 
but only unbounded control are investigated^^ . In other 
words, a comprehensive study of the residual decoherence 
of a system, which is coupled to a general environment 
and to which a general dynamical decoupling is applied, 
is still missing. This is mostly due to the lack of a con- 
venient tool. An unified framework to investigate these 
effects is thus much needed. 

In this work we develop an unified numerical frame- 
work to simulate the renormalization of non-Markovian 



dynamics of a system to which decoupling pulses are ap- 
plied. It is capable of handling both unbounded and 
bounded pulses. It is also able to capture the effect of 
a general environment. Such a framework can help us 
to study the efficiency of a prescribed decoupling pulse 
in realistic situations. It might also open new avenues 
to optimize the decoupling pulses. To see what are the 
necessary ingredients of such a framework, recall that a 
quantum system without decoherence obeys Schrodinger 
equation, resulting a coherent, hence non-Markovian, dy- 
namics. On the other hand, long time dynamics of an 
open quantum system in high temperature regime is typ- 
ically Markovian 15]. A successful decoherence control 
should renormalize the dynamics, moving the system dy- 
namics away from Markovian regime and move it into 
the non-Markovian regime. The residual decoherence, 
however, prevent us from simply using the Schrodinger 
equation to describe the renormalized dynamics. It is 
thus essential that such a framework can describe both 
the Markovian and non-Markovian dynamics for an open 
quantum system coupled to a general environment. Also, 
for numerical purpose a nonconvoluted master equation 
is highly desirable comparing to an integral-differential 
equation. 

The nonconvoluted non-Markovian master equation 
|l6l llTj l derived from the non-Markovian quantum tra- 
jectories [T§ | formalism naturally provides such a frame- 
work. Markovian dynamics and Schrodinger dynamics 
can be seen as two limiting cases of the formalism. The 
formalism has been used to study the decoherence of 
the open quantum system in the non-Markovian regime 
ham as well as the non-Markovian dynamics of a two- 
level atom immersed in a photonic band-gap material p^. 
In some cases the stochastic non-Markovian quantum 
trajectories can be formally averaged, resulting an ex- 
act non-Markovian master equation. In general, when 
the decoupling are applied to the system, an exact non- 
Markovian master equation cannot be formed 16J. But 
an approximated master equation can always be de- 
rived in the nearly Markovian limit |l6j| or weak coupling 
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limit |T3|- It has been shown that the resulting non- 
Markovian master equation is valid in a wide range of 
parameter regime [T^ ■ It is this nonconvoluted master 
equation that will be taken as the starting point of our in- 
vestigation of the renormalization of the non-Markovian 
dynamics. 

The organization of the manuscript is as follows. In 
Section^ we briefly review the non-Markovian quantum 
trajectory formalism and the resulting non-Markovian 
master equation. In section 11111 we derive the non- 
Markovian master equations when the decoupling pulses 
are applied. In Section IIVI we simulate the renormal- 
ized dynamics after decoupling of a system which is 
driven by the Ornstein-Uhlenbech noise while in Sec- 
tion we study the renormalized dynamics of the spin- 
boson model. Three representative spectral densities are 
used. Bang-bang control pulses and continuous decou- 
pling pulses are used in each case. The conclusion in 
Section IVTI brieflv summarizes our findings. 



II. NON-MARKOVIAN QUANTUM 
TRAJECTORIES AND MASTER EQUATIONS 

The total Hamiltonian of our model of the open quan- 
tum system takes the form Htot — Hsys + H^nv + 
Hint 1 where Hsys — Hq + Hc(t) represents the sys- 
tem Hamiltonian which has been split into the orig- 
inal Hamiltonian Hq and the dynamical decoupling 
control Hamiltonian He- Henv = J2x^^^\^>^ repre- 
sents the bosonic Hamiltonian of the environment and 
Hint = (j2x Lg'^a\ + g\a)}j represents the interaction 

between system and environment. In our simulations we 
have set Hq = since we mainly focus on quantum mem- 
ory In this work. 

It has been shown that one can derive linear stochastic 
equations governing the dynamics of the system under 
the influence of the environment Ts'l. The stochastic 
Schrodinger equation can be obtained by expressing the 
total state in terms of the Bargmann coherent state 

of the environment 
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Using this representation the reduced density matrix of 
the system can be expressed as 
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If one choose to use Monte-Carlo integration to evaluate 
this expression, taking over the Gaussian distribution 
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of z— vectors, then r]*{t) can be interpreted as stochastic 
variables with following statistical properties 
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Under this interpretation, Eq.(|21) becomes a stochastic 
Schrodinger equation on the system. However it is still 
difficult to treat Eq.Q due the the functional deriva- 
tive with respect to the noise 7?*(s) under the memory 
integral. A nonconvoluted stochastic equation can be 
obtained by proposing the Ansatz 
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Oit,s,rj{t))\iP{z\t)), 
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where 0{t, s,r]{t)) is a linear operator that has to be 
determined for each case. With this replacement the 
stochastic equation becomes 

j^\ijiz*,t))=[-iHsys+L,^*it)-L^d{t,ij*)] mz*,t)), 

(9) 

where 0{t,r]*) — dsa{t, s)0{t, s,ri*). For some cases 
it is possible to construct exactly the O operator. When 
it is difficult to identify the O operator, it is still possible 
to derive an approximate O operator by perturbation in 
terms of the coupling parameter (171 or the environmen- 
tal correlation time|16i|. If the exact or approximate O 
operator is independent of the noise rj* , a master equa- 
tion can be easily derived and takes the form: 



dpjt) 
dt 



^t[Hsys,p] + [L, p{t)dHt)] + [0{t)p{t), Lt]. (10) 



Since we are interested in the renormalization of the non- 
Markovian dynamics we take the weak coupling limit in 
stead of the nearly Markovian limit. In this limit 
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Hence 
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In the interaction picture with respect to the system 
Hamiltonian p{t) — {t)p{t)U{t), the master equation 
becomes 



dm 

dt 



L(t),p{t) / dsa*{t,s)L\s) 



dsa{t,s)L{s)p{t),L\t) 



(14) 



where L{t) = W{t)LU{t). Note that when a{t,s) 
S(t — s) the master equation reduces to the standard 
Markovian Lindblad master equation in Heisenberg pic- 
ture, i.e. 



dm 

dt 



^[m,mL\t)] + [L{t)m,LHt)]. (is) 



which contains no decoherence terms. The connec- 
tion to the standard dynamical decoupling can be seen 
more transparently if one assumes that a{t,s) — 9{\t — 
s\/Tc)/Tc and Hq = 0. The decoupling conditions then 
becomes 



/ ^ dsUUs)LA,_ 
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(19) 



III. DYNAMICAL DECOUPLING 



which is exactly the standard bang-bang decoupling 
condition 0. 

In the following we will derive the non-Markovian mas- 
ter equations when continuous or bang-bang decoupling 
is applied to the system. Those master equations will 
be used in next two sections to simulate the renormal- 
ized dynamics. First consider the case where L = a^. 
Without the control Hamiltonian the master equation is 



In the language of standard periodic dynamical de- 
coupling one usually starts from the total Hamiltonian 
Htot = ® B^. A periodic control Hamiltonian 

Hc{t) with period Tc is introduced to decouple the sys- 
tem from the environment. Using Magnus expansion [2lj 
the average Hamiltonian governing the stroboscopic dy- 
namics has the form 
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ds 
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A properly designed control Hamiltonian can generate 
an Ui.{t) which results in a zero Htot- Many meth- 
ods have been developed to design the decoupling pulses 
0,0, 0,0, 00- In this work we focus on investigating the 
renormalized dynamics when a decoupling pulses based 
on those methods are applied to the system. It should be 
noted, however, that it is also possible to design decou- 
pling pulses based on Ea. (|14|l . We sketch the procedure 
as follows. First decomposes L into L = La + iLs where 
both La and Lb are Hermitian. It is then easy to verify 
that the decoupling condition (for the case of quantum 
memory) is 

dsa{t, s)La{s) — Xa, / dsa{t,s)LB{s) = Xb, 
10 Jo 

(17) 

where Xa and Xb are real numbers. When this condition 
is satisfied, the master equation turns into 



dp{t) 



dt 



= [La + iLB, mi^A - iXB)] + h.c. 

= -t[i-2X aL b), m] - i[i+^^B La), m] 
= -z[#e//,p(t)], (18) 
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where we have decomposed a(t, s) into its real part 
/x(t, s) and imaginary part v{t, s) but with the conven- 
tion a = p, — iv. Note that y{t, s) doesn't contribute 
the decoherence in this case. According to Ref.Q, a con- 
tinuous decoupling can be achieved by using the control 
Hamiltonian Hc{t) = ax{t)ax, where we have assumed 
Hq — since we are interested in the quantum memory. 
Under such a decoupling the master equation becomes 



dm 

dt 



(Tz(t),m I dsa*{t,s)az{s) 





h.c. 



(21) 



To proceed, one needs to find the associated time evolu- 
tion operator 



TT /i\ f dua:c(u)(7-j. 

Uc{t) Jo = 

where A{t) — duax{u), and 



cos{A{t)) -ism{A{t)) 
4sm{A{t)) cos{A{t)) 
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cos(2yl(s)) -ism{2A{s)) 
ism{2A{s)) cos(2yl(s)) 



(23) 

It has been shown in Ref . 0] that decoupling condition 
can be achieved if 2A{Tc/2) = tt and 2^(rc/2 + s) = 
TT + 2A{s). Under this condition one has 
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i sin(2A(s)) - i sin(2A(s)) cos(2A(s) - cos(2A(s)) 
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(24) 



Hence at the ideal limit where Tc 
totally suppressed, i.e., 

dm 

mil 

Tc 



the decoherence is 



lim 

*0,N^oo,t=NTc 



dt 



0. 



(25) 



t=NTc 



However when is finite, the decoupling cannot be per- 
fectly achieved. The error associated with finite Tc will 
depend on the correlation function a{t, s) and the enve- 
lope function a^it). The master equation associate with 
a finite T^. can be easily derived. After some algebra we 
find 

^ = 2sin(2^(t)) / dsiJ.{t,s)sm{2A{s)) 
dt Jq 



Pii(i) -Poo(t) 
-Poi{t) - Pio{t) 
ft 



-Poi(i) - Pio(i) 
-Pn{t) + Poo{t) 



+ 



2sin(2^(t)) [ dsu{t, s) cos{2A{s)) 
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Poi(i) + Pio(t) -Poo(i) + Pn{t) 
-Poo(i) + Pn{t) -poi(i) - Pio(i) 



(26) 



where we have used a{t, s) = ii{t, s) — iv{t, s). Once the 
envelop function ax{t) is identified, this equation can be 
used to calculate the renormalized dynamics. Note that 
under continuous decoupling both ^, and v contribute to 
the decoherence. 

On the other hand for bang-bang control the function 
A(t) takes the special form 



1=0 



(27) 



It is then straightforward to show that ^^(s) = /(s)^^ 
where /(s) = +1 for s € (2i^,2i + 1^) and /(s) = 
-1 for s e {2i + 1^, 2i + 2^), i = 1,2,..., [^]. The 
resulting master equation is 



dm 

dt 



dsiJ,{t, s)f{s) 
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In other words, applying bang-bang control results in the 
renormalization of the function fi{t, s). The renormalized 
/i(t, s) can be evaluated numerically as follows 



/ dsii{t,s)f{s) 
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= V / dsM(t,s)(-l)'=-V, 
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(29) 
dsii{t, s)az. 



We observe that the master equation under bang-bang 
control takes a simpler form compared to the one under 
continuous decoupling. This suggests that there might 
be some dynamical effects which cannot be captured if 
one only studies the ideal bang-bang limit. 

Consider next the case where L = a-. Without the 
control Hamiltonian the master equation is 



dm 

dt 
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+ipio 



(30) 



Note that in this case both fi and u contribute to the 
decoherence. When continuous decoupling is applied the 
master equation becomes 



dp 
'dt 



(31) 
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and 



i>(t) = 4e-^2A(t) t dsv{t,sy^^^'\ (33) 
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While when bang-bang control is applied, the master 
equation becomes 



dm 
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In contrast to the case where L = cr^ , in this case both p, 
and V contribute to the decoherence. 



IV. ORNSTEIN-UHLENBECK NOISE 

In this section we study how the dynamics of a system 
driven by the Ornstein-Uhlenbeck noise is renormahzed 
by the bang-bang and continuous decouphng pulses. The 
Ornstein-Uhlenbeck noise is characterized by the expo- 
nential correlation function 



a{t,s) = n{t,s) = ^e" 



(35) 



This particular form of the correlation function enables 
us to carry out analytically most of the calculation. The 
main purpose is to illustrate how the decoherence is 
renormalized by the bang-bang or continuous decoupling 
pulses. First consider the case where L = a^- Using 
the result in proceeding section, we find that without the 
control Hamiltonian the off-diagonal term of the density 
matrix decays as follows: 
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dt 



i = -2 (l - e-*/- 
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When a control Hamiltonian of the form Hc{t) = ax{t)(Tx 
is applied and the resulting A{t) has the form A{t) — 
■nt/Tc, the renormalized master equation becomes 



dp{t) 2Af,(t)sin(27ri/Tc) 



dt 



1 + {2nT/T,y 



(37) 



pii(t) - poo(i) ~PQi{i) - Pw{t) 
-poi{t) - Pio{t) -pii{t) + poo{t) 



where 



M,(t) = e-*/^(^)-^cos(^)+sin(^). (38) where 
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It is important to observe that 

2Ms{t) sin {2TTt/Tc) 
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Ta^O.N^oo,t=NTc 



1 + {2^T/T,y 



= 0. 



(39) 



If a different envelope function a^it) is used, the resulting 
renormalized master equation will have a similar form, 
but a different time-dependent pre-factor will appear. 

On the other hand, if bang-bang decoupling is applied 
the renormalized master equation becomes 



dpai 



= -2e- 



dt 

Note again that 
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2e~ — 
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0. 



V l + e^ / 

(41) 

Consider next the case where L = (7+ . In this case the 
master equation without the control is 



2/5ii(t) -2poi 
-2pQi -2pn{t) 



(42) 
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FIG. 1: Time evolution of poi(i) for a system driven by the 
Ornstein-Uhlenbeck noise and bang-bang or continuous con- 
trol. In all figures we set r — 0.5. Upper-left: L — az and 
bang-bang control is applied. Bottom-left: L — and con- 
tinuous control is applied. Right: L = and continuous 
control is applied. In this case poi(i) is complex and its ab- 
solute value and real part are plotted respectively. 



When the control Hamiltonian Hc{t) = az{t)az is applied 
and A{t) — nt/Tc, the renormalized master equation is 



dm 



dt 



(43) 



-(M(0 + Ai*(i)Poo(i) -M*(OPoi(i) 
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While when bang-bang control is applied, the resulting 
master equation is 



dm 

dt 
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„ Ta 1 — e' 

1 -I- 2e— 



2pii(i) -Poi(i) 
-pio{t) -2pn{t) 



(45) 



It is instructive to simulate numerically the time evo- 
lution of PQi{t) under bang-bang and continuous decou- 
pling. We assume that initially the system is at pure state 
IV^) = ^(|0) + |1)). In FigUKleft) we plot pai{t) as func- 
tion of time for the case where L = Gz- The renormal- 
ized dynamics under bang-bang control and continuous 
decoupling are plotted respectively. In this case PQi{t) is 
always real. In Fig^ (right) we plot the real part and the 
absolute value of poi {t) for the case where L — a-. We 
only plot the renormalized dynamics under continuous 
decoupling since the dynamics under bang-bang control 
is identical to the case where L = az with a different 
overall constant. 
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FIG. 2: fif '^{t) as a function of time with an Ohmic envi- 
ronment and A — 20. (a) Without bang-bang control, (b) 
Tc = 0.5. (c) Tc = 0.25. (d) T,. = 0.125. (e) n = 0.0625. 



FIG. 3: Ui '^{t) as a function of time with an Ohmic envi- 
ronment and A = 20. (a) Without bang-bang control, (b) 
Tc = 0.5 (c) Tc = 0.25. (d) Tc = 0.125. (e) Tc = 0.0625. 



It is evident that when bang-bang or continuous decou- 
pling is applied the decoherence is suppressed and better 
result is achieved for shorter Tc. It is important to note 
that the efficiency for those two class of decoupling is dif- 
ferent. Continuous decoupling results in a much better 
suppression. The continuous decoupling can already sup- 
press the decay of poi when Tc r. Although we only 
show one particular realization of the continuous pulse, 
numerical evidence shows that continuous pulse always 
results in better suppression for this noise model. Fur- 
ther optimization on pulse shape is possible and worth 
more investigation. It should also be noted that that (Tz 
and coupling to the environment results in quantita- 
tively different renormalized dynamics. 



A. Ohmic and Supra-Ohmic 

Consider first the Ohmic and supra-Ohmic environ- 
ment. It is instructive to exam how and v are renor- 
malized by the bang-bang control. In Fig[21 we plot 
11^'^ (t) = rfs/i(t, s)/(s) for an Ohmic environment. 
It is clear from the figure that when Tc is not small 
enough, fi^'^'{t) might be larger than the original fi{t) 
in some temporal regime, which is undesirable. For small 
enough Tc the ii{t) is always suppressed. In FigElwe plot 
j^B.B. ^ J* (^gp^f^ s)/(s) for the same system. We find 
that lyit) is always suppressed by the bang-bang control. 
Continuous decoupling, on the other hand, will mix the 
real and imaginary part of a{t, s). In some cases this will 
deteriorate the efficiency of the decoupling. In the follow- 
ing we present the numerical results for various scenarios. 



SPIN-BOSON MODEL 



1. L = az, bang-bang decoupling 



In this section we turn our attention to the spin-boson 
model. It is known that the zero temperature correlation 
function of bosonic environment can be written as 



a{t, s) = / dujJ{Lu)[cos{Lu{t — s))—ism(Lij{t~s))], 
Jo 

(46) 

where the spectral density J{u!) has the form 



J(w) = w^e"'^/^^^. 



(47) 



When p = 1 it corresponds to the Ohmic environment, 
when p > 1 it corresponds to the supra-Ohmic environ- 
ment, and when p < 1 it corresponds to the sub-Ohmic 
environment. For 1// type sub-Ohmic environment, i.e., 
p = —1, an infrared cutoff A/^ is also necessary. 



When L = (7z the master equation at zero temperature 
without the control Hamiltonian is reduced to 

= -4 f ds^i^,^{t - s)poi(i), (48) 



dt 



where 



At^s)= / duujPe-'^''^cos{uj{t~ s)). (49) 



When bang-bang control is applied the master equation 
reads 

= -4^* dsMp,A(i - s)/(s)poi(i), (50) 

We define T2(A) to be the time which satisfies the con- 
dition /5oi(T2) = e~^poi(0)- In the Markovian limit this 
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FIG. 4: poi(i) under bang-bang control for a system coupled 
to an Ohmic (p = 1) and Supra-Ohmic (p = 3) environment 
and L — (Jz- Solid line represents the free decay without 
bang-bang control. 



FIG. 5; poi{t) under continuous decoupling for a system cou- 
pled to an Ohmic (p = 1) and Supar-Ohmic (p = 3) environ- 
ment and L = CTz . Solid line represents the free decay without 
decoupling control. 



definition will coincide with the typical definition of T2. 
In Fig0]we plot the time evolution of poi{t) under bang- 
bang control for a system coupled to an Ohmic {p = 1) 
and a supra-Ohmic {p ~ 3) environment. We have picked 
two different cut-offs A for each environment. It is evi- 
dent that the free decay is not a simple exponential and 
the detail functional form of the decay depends on p and 
A. Since the decay is non-exponential, there is no sin- 
gle time scale to which Tc can be compared. Roughly 
speaking, however, one can still say that the decoupling 
is efficient when Tc is smaller than T2(p, A). 



3. L = (J-, bang-hang decoupling 

Consider next the case where L = a^. Without the 
decoupling pulse the differential equation for poi(i) reads 

^^d/ ^ ^ J ^'^ 
while when bang-bang control is applied, it becomes 



2. L = Gz, continuous decoupling 

When continuous decoupling is applied the differential 
equation for pai{t) reads. 



dpoiit) „ . ,2TTt. /■* , , . ,2ns. 

— ^ = -2sm( — ) / dspp^A[t - s)sm( — )poi, 

(51) 

where we have assumed the same initial condition and 
A{t) — nt/Tc- Note that in this case poi{t) remains real 
during the evolution. In Fig|31 we plot the time evolu- 
tion of poi{t) under continuous decouphng control for a 
system coupled to an Ohmic {p = 1) and a supra-Ohmic 
{p = 3) environment. Similar to the case where the sys- 
tem is driven by the Ornstein-Uhlenbeck noise, continu- 
ous decoupling is more efhcient for the same Tc- 



-^^^-f^ ^ ^ J '^^ (Aip,A(i, s) + iVp.Ait, s)) f{s)poi{t), 

(53) 

In contrast to the case where L = CTz , it is now necessary 
to evaluate the renormalization of both p and v. In FigEl 
(Figl?!) we plot the renormahzed dynamics of poi{t) un- 
der bang-bang control for a system coupled to an Ohmic 
(supra-Ohmic) environment. The behavior of the real 
part of poi (t) is qualitatively similarly to the case where 
L = Uz- The imaginary part, on the other hand, deviates 
from zero as time goes. In most cases it grows monoton- 
ically. In some cases, however, it shows quasi-oscillation 
behavior when no decoupling pulse is applied. For some 
range of the Tc, the bang-bang control cannot efficiently 
suppress the growth of the imaginary part. For small 
enough Tc, the renormalized dynamics of the imaginary 
part resumes its monotonic increase, with a suppressed 
increase rate. It should be noted that the overall effi- 
ciency is qualitatively similar to the case where L = Uz 
since the decay of the real part of poi dominates the de- 
coherence. 
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FIG. 6: poi{t) under bang-bang control for a system coupled 
to an Ohmic {p — 1) environment and L = a-^. Solid line 
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FIG. 7: poi(i) under bang-bang control for a system coupled 
to a Supar-Ohmic {p — 3) environment and L = a-. Solid 
line represents the free decay without bang-bang control. 
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FIG. 8: poi{t) under continuous control for a system coupled 
to an Ohmic (p = 1) environment and L = a-. Solid line 
reoresents the free decav without bane-bane control. 
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FIG. 9: poiit) under continuous control for a system coupled 
to a Supar-Ohmic {p — 3) environment and L = a-. Solid 
line represents the free decay without bang-bang control. 



4-. L = 0-, continuous decoupling 

When continuous decoupling is applied the differential 
equation for PQi{t) reads, 

(54) 

In Fig|Sl (FigO we plot the renormalized dynamics of 
Pmit) under continuous control for a system coupled to 
an Ohmic (supra-Ohmic) environment. The salient fea- 
ture to be observed is the dynamics of the imaginary part 
of poi(*)- It only increases slowly for the case of free de- 
cay. When the continuous decoupling is turned on, how- 
ever, the imaginary part grows much more rapidly. This 
behavior results in the deterioration of the efficiency. The 
growth of the imaginary part can be suppressed by de- 



creasing the Tc. But a very small is needed to suppress 
the imaginary part to the same level as the one in free 
decay. This feature is due to the fact that when continu- 
ous decoupling is turned on, /i and v are not separately 
renormalized. The decoupling pulse will mix /i and 
resulting a more rapid growth of the imaginary part. 



B. Sub-Ohmic 

Now turn our attention to the sub-Ohmic environment. 
We will focus on the 1// type environment. For 1// type 
sub-Ohmic environment the correlation function takes 
the form 
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FIG. 10: poi{t) under bang-bang (upper half) or continuous 
(bottom half) control for a system coupled to an 1//- noise 
environment and L — az- Solid line represents the free decay 
without decoupling control. 



FIG. 11: poiit) under bang-bang (right two) or continuous 
(left two) control for a system coupled to an l//-noise en- 
vironment and L = <j-. Solid line represents the free decay 
without decoupling control. 



a[t,s)= I duj [cos{uj{t ~ s)) — ism{uo{t — s))] ^ 

Jair ^ 

(55) 

where an infrared cutoff is introduced to ensure the con- 
vergence. In Fig llOl we plot the renormaUzed dynamics 
of pqi (t) for a system with L = az under bang-bang or 
continuous decoupling. We observe that continuous de- 
coupling is more efficient. In Fig llll we plot the renor- 
malized dynamics of poi(^) for a system with L = u- 
under bang-bang or continuous decoupling. In this case 
we observe that imaginary begins to grow rapidly once 
the continuous decoupling is turned on. This is again 
due to the mixing of //(t) and ^{t) via the continuous de- 
coupling. In this scenario the bang-bang control is more 
efficient. 



VI. SUMMARY AND DISCUSSION 

In summary we have developed a numerical framework 
to investigate the renormalization of the non-Markovian 
dynamics due to the dynamical decoupling pulses. The 
nonconvoluted non-Markovian master equation used in 
this work is derived from the non-Markovian quantum 
trajectories formalism. The resulting master equation is 
written in a convenient form such that once the decou- 
pling pulses are prescribed, the renormalized dynamic 
can be readily simulated. In order to validate the frame- 
work we have performed a comprehensive simulation of 
the renormalized dynamics under following scenarios. We 
have investigated two different system-environment inter- 
actions, namely L — Oz and L = a^. Three representa- 
tive environmental spectral densities are used, including 
supra-Ohmic, Ohmic, and sub-Ohmic environment. Two 



kinds of decoupling pulses are considered. One is the 
typical bang-bang decoupling while the other one is the 
continuous decoupling. Note also that the period of the 
decoupling cycle, T^, is finite in all of the simulations. 

The numerical results presented in Section IIVI and 
clearly demonstrate that the framework is capable of sim- 
ulate the renormalization of the dynamics of a system 
coupled to a general environment under the action of 
a general decoupling pulse. The renormalized dynam- 
ics, in turn, determined the true efficiency of a decou- 
pling pulse. We observe that different decoupling pulses 
sometimes give rise to qualitatively different renormal- 
ized dynamics, and the difference depends on both the 
system-environment interaction and the environmental 
correlation function. Recall that, however, these decou- 
pling pulses are designed solely based on the form of the 
system-environment interaction, and they all attain the 
desired decoupling in the ideal limit. To better under- 
stand the implication of this observation one should no- 
tice that bang-bang decoupling can be viewed as an un- 
physical limit of the continuous decoupling Q. Within 
the framework presented in this work, bang-bang decou- 
pling simply corresponds to a very special envelop func- 
tion. On the other hand, only one particular envelop 
function A{t) is used for continuous decoupling in this 
work. In other words, two special envelope function are 
chosen to carry out the simulations throughout this work. 
In reality the real envelop function of the pulse, which 
might be designed based on other methods, should fall 
somewhere in between these two limiting cases. The pro- 
found difference in the renormalized dynamics observed 
in this work thus indicates that the real efficiency de- 
pends crucially on the actual shape of the envelop func- 
tion, the spectral density of the environment, and the na- 
ture of the system-environment interaction. For example, 
our results show that continuous decoupling is more ad- 
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vantageous in many cases, as can be seem in Fig^ Fig01 
Fig|Sl and Fig^| However when L = tT_ the continu- 
ous decoupHng will induce mixing between the real and 
imaginary part of the environment correlation function, 
making the decoupling less efficient. Above observation 
actually leads to the possibility of further optimizing the 
decoupling pulses basing on the master equation Eq H14|) 
and using techniques such as pulse shaping 22]. 

The framework developed here can server multiple pur- 
poses. The primary goal is to server as a convenient tool 
in investigating the efficiency of a prescribed decoupling 
pulse. Although we have only shown the cases where the 
deterministic, periodic pulses are applied. The frame- 
work can also easily simulate the reriormalized dynamics 
induced by a random decoupling 5, 6] or a concatenated 
dynamical decoupling ^J. It can also be used as an 
alternative starting point to design the decoupling pulse. 
One can naturally view the non-Markovian master equa- 
tion under decoupling as a dynamical map, whose con- 
vergence the ideal fixed point of zero system-environment 
interaction can be studied. The framework is not limited 



to a single two-level system coupled to a general environ- 
ment. It is also suitable for studying the renormalized dy- 
namics of a multi-level system. Such a multi-level system 
can represent a multi physical or logical qubit system. 
When necessary, one can unravel the master equation and 
use the corresponding stochastic Schrodinger equation, 
which is numerically advantageous for a large system. It 
is also important to note that if eventually a continuous 
measurement interpretation of non-Markovian stochastic 
equations is developed, this framework can be adapted 
to study the error correction and feedback control in the 
non-Markovian regime. 
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